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ABSTRACT 

We present the catalog of the Swift X-ray Cluster Survey (SWXCS) obtained us¬ 
ing the archival data of the X-ray Telescope ( XRT) onboard the Swift satellite acquired 
from February 2005 to November 2012, extending the first release of the SWXCS. The 
catalog provides positions, soft fluxes and, when possible, optical counterparts for a 
flux-limited sample of X-ray group and cluster candidates. We consider the fields 
with Galactic latitude \b\ > 20° to avoid high HI column densities. We discard all 
the observations targeted at groups or clusters of galaxies, as well as particular ex- 
tragalactic fields not suitable to search for faint extended sources. We finally select 
~ 3000 useful fields covering a total solid angle of ~ 400 deg 2 . We identify extended 
source candidates in the soft-band (0.5-2 keY) images of these fields, using the soft¬ 
ware EXSdetect, which is specifically calibrated on XRT data. Extensive simulations 
are used to evaluate contamination and completeness as a function of the source sig¬ 
nal, allowing us to minimize the number of spurious detections and to robustly assess 
the selection function. Our catalog includes 263 candidate galaxy clusters and groups, 
down to a flux limit of 7 x 10 -15 erg cm” 2 s _1 in the soft band, and the logN-logS is 
in very good agreement with previous deep X-ray surveys. The final list of sources 


'CAS Key Laboratory for Research in Galaxies and Cosmology, Department of Astronomy, University of Science 
and Technology of China, 230026, Hefei, Anhui, RR. China; liuteng@ustc.edu.cn 

2 INAF, Osservatorio Astrofisico di Firenze, Largo Enrico Fermi 5,1-50125, Firenze, Italy 

3 INAF, Osservatorio Astronomico di Brera, Via Brera 28,1-20121 Milano, Italy 

4 Universita degli Studi di Ferrara, Dipartimento di Fisica e Scienze della Terra, Via Saragat 11-44121 Ferrara, Italy 

5 INAF, Osservatorio Astronomico di Brera, Via Bianchi 46,1-23807, Merate (LC), Italy 

6 University of Massachusetts, Department of Astronomy, LGRT-B 619E, 710 North Pleasant Street, Amherst, MA 
(USA) 




- 2 - 


is cross-correlated with published optical, X-ray, and Sunyaev-Zeldovich catalogs of 
clusters. We find that 137 sources have been previously identified as clusters in the 
literature in independent surveys, while 126 are new detections. At present, we have 
collected redshift information for 158 sources (60% of the entire sample). Once the 
optical follow-up and the X-ray spectral analysis of the sources are completed, the 
SWXCS will provide a large and well-defined catalog of groups and clusters of galax¬ 
ies to perform statistical studies of cluster properties and tests of cosmological models. 

Subject headings: galaxies: clusters: general - cosmology: observations - X-ray: 
galaxies: clusters - surveys - catalogs 


1. INTRODUCTION 


Groups and clusters of galaxies are the most massive, gravitationally bound structures in 
the Universe and their hot Intra Cluster Medium (ICM) make them to appear as prominent ex¬ 
tended sources in the X-ray sky. Therefore, X-ray cluster surveys are among the most efficient 
tools to constrain cosmological parameters and primordial density fluctuations. A large and com¬ 
plete catalog of groups and galaxy clusters spanning a wide range of redshifts would be crucial 
to make s ignificant steps forward towards the understanding of cosmic structure formation and 
evolution (RosatietaL 2002a : Schuecker 2005 : Voit 2005 : Borgani 2008 ), t he ch e mical and ther 


2008 


(see 


modynamical cosmic histor y of the ICM (lEttori et al.ll2004l: iBalestra et al.l 120071: iMaughan et al. 


Anderson et al. 


Vikhlinin et al. 


2009i 


2009 ), and to provide an accurate measu rement of cosmological parameters 


Mantz et al.ll2010l: I Allen et al.ll2011 ). 


However, this task is not within reach of current X-ray missions. All the major X-ray facilities 
existing today have not been designed for surveys, and have a low efficiency in discovering rare 
objects like galaxy clusters, particularly at high redshifts. The main characteristics required for an 
effective X-ray survey mission for extended sources are: large field of view (FOV, of the order of 
1 deg 2 ), high angular resolution (of the order of few arcsec), low background, and a large effec¬ 
tive area (of the order of 10 4 cm 2 ). Loo king at the near future, the upcoming mission eROSITA 
(IPredehl et al.ll2010l: iMerloni et al.l 120 1 2li will finall y provide an X-ray all-sky coverage 20 years 
after the ROSAT All Sky Survey (lYoges et al .111 99% . down to limiting fluxes more than one order 
of magnitude deeper than ROSAT for extended sources. Therefore eROSITA will considerably 
increase the number of X-ray groups and clusters particularly at low and moderate redshifts. How¬ 
ever, its limiting flux is predicted to be ~ 3.4 x 10" 14 erg cm" 2 s -1 after four years of operation, well 
above the level below which the majority of the high-z clusters, and medium and high-z groups, are 
currently found. In addition, its low effective area above 2 keV severely limits the characterization 
of the ICM in high temperature clusters (see|Borm et al. 2014). 
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At present, the best resource to build X-ray cluster samples is provided by the still increasing 
archives of the major X-ray facilities, Chandra and XMM-Newton. For a revie w of the ongoing X- 
ray cluster surveys, updated to the year 2012, see Table 1 in Tun do et al. (2012). In this framework, 
we recently presented the Swift X-ray Cluster Survey (SWXCS iTundo etah 2012 . herea fter Paper 
I), which is based on the archival data of the X-ray telescope (XRT. iBurrows et al.ll2005l) onboard 
the Swift satellite ( Gehrels et a l. 20040 . In spite of its small collecting area (about 1/5 of Chandra ) 
XRT has two characteristics which make it an efficient instrument for detection and characteriza¬ 


tion of extended sources: a low background (IMoretti et al. 2009 ) and a constant a ngular resolution 
(with a Half Energy Width HEW = 18") across the entire FOV ( Moretti et al. 2007 ). We note 
that XRT angular resolution is as good as the resolution of XMM-Newton at the aimpoint, and 
therefore better than XMM-Newton when averaged over the FOV. The first catalog of the SWXCS 
project, including 72 clusters and groups, has been presented in Paper I, while the X-ray spectral 
analysis for more than half of this sample is presented in Tozzi et al. (1201 4 hereafter Paper II). 


In Paper I we used only the Gamma Ray Burst (GRB) follow-up observations of XRT released 
before April 2010. The sample we built in Paper I has shown the efficiency of an X-ray telescope 
as small as XRT in finding and characterizing X-ray extended sources. The natural next step is 
the inclusion of the entire Swift-XRT archive, which is the goal of this work. With this aim, we 
developed a software designed for the detection and photometry of extended sources and optimized 
for the characteristics of XRT data ( EXSdetect , Fiu et al.ll2013 ). The source detection method used 
in EXSd etect, is a combination of Voronoi Tessellation and friend-of-friend algorithms (VT+FOF, 
Ebeling & Wiedenmann 1993 ). This method does not require a priori assumptions about the shape 
and size of the sources, and it is particularly efficient when applied to X-ray images which are 
characterized by many empty pixels. Its main limitation consists in the blending effect of merging 
neighboring sources into one. Spurious extended sources may occur due to the bridging of two or 
more faint unresolved sources, or to the overlap of true extended emission with the wings of bright, 
unresolved sources. To mitigate this effect, in EXSdetect we developed an accurate deblending 
procedure which is very effective in identifying and separating most of the blended sources, and 
eventually removing the unresolved sources mistakenly included within extended emission. 


The efficiency of our detection algorithm as a function of the exposure time and the source 
flux, and the accuracy in the source photometry, are investigated by extensive imaging simulations. 
Most importantly, our simulations allow us to evaluate the contamination (number of spurious 
extended sources) and the completeness of the catalog as a function of the source flux. In the 
simulations we make use of an empirical model of the point spread function (PSF) of XRT. Thanks 
to our simulations, we can identify an optimal threshold in the source photometry above which our 
catalog reaches the required completeness and purity. This threshold directly provides a position- 
dependent flux limit for each field, and hence a selection function depending on the physical source 
flux for the entire survey. This step is particularly relevant since, once the selection function and 
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the contamination level are accurately predicted, the sample can be used for statistical studies. 

This work is the extension of the previous SWXCS catalog (Paper I) to the entire Swift-XRT 
archive as of November 2012, which includes more than 10000 fields, as opposed to the ~ 300 
fields used in Paper I. In addition, we apply for the first time the EXSdetect software to the XRT 
data to achieve better accuracy and sensitivity. The ultimate goal of SWXCS is to provide a well- 
defined, large catalog of X-ray selected groups and clusters to investigate X-ray properties and 
perform classical cosmological tests. The Paper is organized as follows. In we describe the 
selection of the Swift-XRT fields suitable for our serendipitous survey, and describe the identifica¬ 
tion and classification of the extended sources. In CU we present the SWXCS catalog and compare 
it with previous works in the field. In §|4] we compare this catalog to the previous release. In §[5] 
we search for counterparts in the optical, X-ray, and Sunyaev-Zeldovich (SZ) catalogs for all our 
sources. Finally, our findings are summarized in ® 


2. FIELD SELECTION AND SOURCE IDENTIFICATION 
2.1. Field Selection 


From the entire Swift XRT archive in the period February 2005 - November 2012, we select 
all the fields which can be used to build an unbiased, serendipitous X-ray cluster catalog. Firstly, 
we exclude all the fields whose Galactic latitude \b\ < 20°. Although these fields could in principle 
be used to search for bright extended sources, they are typically very crowded, which would cause 
severe blending problems when spatial resolution is limited. Moreover, a significant fraction of the 
soft band emission from groups and clusters would be absorbed by the high HI Galactic column 
density. The search for the brightest groups and clusters in the Galactic fields will be performed 
with a different technique in a dedicated paper (Moretti et al. in preparation). 


Secondly, we exclude the shallow fields whose exposure time is t exp < 3000 s. This limit 
guarantees > 100 total photons in the soft band in each field (see Figure d). This threshold rep¬ 
resents the minimum number of photons to sample the background in an X-ray image, which is 
a critical step to identify extended source as enhancement of the photon density with respect to 
the background level. We set a lower limit of 100 to the number of total photons based on exten¬ 
sive simulations of background-only ima ges, wher e we tested the capability of recovering the true 
background in our algorithm (see Liu et al. I2013L and Appendix A). This forces us to discard a 
large number of fields (see lines in Figured), which, however, would have contributed only at very 
high fluxes (> 10 -12 erg cm" 2 s _1 ), where the expected number of clusters is very low. 


We are also forced to exclude the deepest field (164440+573434) whose exposure is 1.9 Ms. 
This choice, which appears to be particularly unfortunate, is due to the large background in the fi- 
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nal image. While the internal parameters of EXSdetect are optimized for almost the entire range of 
exposure among our fields, EXSdetect becomes unstable when applied to XRT images with expo¬ 
sure times around 2 Ms. Therefore, we choose to discard the field 164440+573434, which would 
require a different approach with respect to the rest of the survey. We note that the performance 
of EXSdetect in the second deepest field, whose exposure time is 1.1 Ms, has been successfully 
tested with our simulations. 



Exposure time (s) 


Fig. 1.— Number of filled pixels (N) in the soft-band image versus exposure time (/) for each field 
of the XRT archive as of December 2012. The red line (N = 0.03/) shows the lower envelope of 
the distribution. The horizontal blue line shows the condition N > 100 which ensure a reliable 
background measurement, while the vertical black line shows the conservative threshold t exp > 
3000 which guarantees the condition N > 100. Note that the large majority of fields with t exp < 
3000 and N » 100 are dominated by one or few bright sources, therefore, despite the large 
number of filled pixel, their background is poorly characterized. 

A further, fundamental step is to filter out all the fields which are directly targeting groups 
or clusters of galaxies, since they would clearly introduce a positive bias towards the detection of 
extended sources. In principle, one can simply excise the targeted group or cluster, and use the 
rest of the field. However, due to the correlation function of dark matter halos, to have an unbiased 
sample one should also exclude extended sources with similar redshifts. Since this information 
is not available for many of our sources, we decide to remove all the fields which, on the basis 
of the target name and coordinates, are aiming at groups and clusters of galaxies. We also filter 
out observations targeting nearby galaxies, because such galaxies appear as bright extended X-ray 
sources, whose emission is not associated to ICM but mostly to X-ray binaries and massive star 
formation events. On the other hand, all the fields targeting at GRBs and AGNs are included in 
this survey. In Paper I we already showed that GRBs show no spatial correlation with galaxy 
clusters, neither are AGNs expected to be correlated with clusters. It is actually expected that 
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AGNs are suppressed in cluster environments, at least locally (Khabiboulline et al, 2014). Other 
multiwav elength studies on th e occurrence of AGNs in clusters provide inconclusive results (for 
exam ple, Pimbblet et al.112013l : iKoulouridis et al.112014 Ehlert et al.112013c iKlesman & Sarajedini 


2014). Recently, it has been found that the fraction of luminous AGNs in clusters reaches that 


in the field in the redshift range 1 < z < 1.5 (Martini et al . 2013) . Therefore, we conclude that 
no significant bias in cluster detection is expected from the inclusion of all the fields targeting at 
AGNs. 


The task described above is not straightforward, since the target information in the header of 
Swift-XRT event files is often different from the standard naming conventions or incomplete. We 
go through the keywords of target names and coordinates in all the fields, and identify the targets 
in the NASA Extragalactic Database (NED) when needed. Here we describe the selection rules 
applied to the Swift-XRT archive in order to identify an unbiased subset of fields in details. Only 
fields which survived the first triage, i.e., with Galactic latitude \b\ < 20° and t exp > 3000 s, are 
considered here. For completeness we also list the excluded fields. 


Selected fields: 


486 GRB follow-up fields, including all the fields previously used in Paper I; 

698 fields whose targets are found in the NED and are classified as AGN; 

654 fields whose targets match an AGN within 5"of their coordinates in the NED; 

22 Swift-BAT triggered observations, whose targets are variable hard-X-ray sources; 

136 fields targeting Fermi/LAT gamma-ray sources (c orre sponding to header keywords: 0FGL, 
1FGL and 2FGL lAbdo et al.ll201ol: iNolan et alfioTj) ; 


71 fields targeting INTEGRAL gamma-ray sources (header keyword IGR. Bi rd e t al. 20041) ; 

84 fields targeting ROSAT detected sources which are not classified as galaxy clusters in NED; 

401 fields whose target names are found in the NED, and are not classified as galaxy clusters or 
groups; 


422 fields targeting Swift-BAT detected sources (mostly local AGN); 

10 safe pointings, which are carried out when the telescope looses coordinates; 
12 fields targeting pulsars. 
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Excluded fields: 


55 fields targeting Abell clusters, Hydra cluster, Coma and the Crab Nebula; 

27 fields targeting the following nearby galaxies: M31, M33, M63, M67, M81, M82, M87, M100; 
17 fields targeting comets, which may show diffuse X-ray emission; 

157 fields targeting nearby galaxies cataloged as NGC or Mrk; 

142 fields within ll°of the Large Magellanic Cloud (LMC) or 6°of the Small Magellanic Cloud 
(SMC); 

111 fields targeting supernovae, which are often hosted by nearby galaxies; 


26 fields targeting Hickson Compact Groups (HCG, Hickson 119821) : 

2 fields targeting ROSAT detected sources which are classified as galaxy clusters in the NED; 
30 fields whose target names are found in NED, and are classified as galaxy clusters or groups; 


9 fields targeting Redshift Survey Compact Groups (RSCG, Barton et al. 1996 1: 

52 fields targeting other sources related to clusters in dedicated observational programs; 
8 fields significantly affected by stray-light; 

762 fields with unknown target classification. 


As previously mentioned, we still have a large number of fields (more than 700) targeting at 
sources whose classification is uncertain. We conservatively discard all these fields with unknown 
classification, although many useful fields may be lost with this choice. However, this enables us to 
avoid any field selection bias, which is a critical requirement for statistical studies and cosmological 
tests. For consistencies with the first SWXCS catalog, we include 8 fields which would be excluded 
with these selection rules (because inside the LMC region) but have been used in Paper I. This 
inclusion does not have any significant effect on the final catalog. Finally, we select 3004 fields 
which provide a truly serendipitous sampling of the extragalactic sky. The positions in the sky of 
the aimpoints of the selected fields are shown in Figure O 
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Fig. 2.— Positions of the 3004 selected fields of the SWXCS. Red points show the GRB fields 
used in Paper I, while blue points show the new fields used in this work. The positions of LMC 
and SMC are marked with large dots. 


2.2. Selection Of Group And Cluster Candidates 


The XRT data reduction is described in Paper I. We consider only the soft-band (0.5-2 keV) 
images for source detection, since the inclusion of the Swift-XRT hard-band images is not useful 
to identify nor to characterize the detected sources. 


At variance with Paper I, where the sources were identified on the basis of a growth-curve 
method, we use here for the first time the EXSdetect so ftware. This algorithm and its performance 
on the Swift-XRT data are described in full details in Liu et al. (120131) . However, since the soft¬ 
ware, which is publicly available and currently updated on the SWXCS website (http: //www. arcetri . astro. i 
and http://swxcs .ustc. edu.cn), evolved significantly in the meanwhile, in the Appendix we 
describe in detail the most relevant changes in the current version (v3.0). The first important 
change concerns the evaluation of the background (see Appendix lAl). while the second concerns 
the source classification scheme (see Appendix©. Another modification introduced in this ver¬ 
sion, is a different treatment of the sources at large off-axis angles (6 > 9'). This is necessary since 
the in-flight calibration of the Sw ift-XRT PSF shows that at 6 ~ 10' the HEW increases signifi¬ 
cantly by ~ 40% (Morctti et a l. 20051) . Therefore we run EXSdetect with a different PSF model to 
match the expected behavior of the PSF at large off-axis angles, only for the source candidates at 
d > 9'. 
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We achieve a strong control on the purity and the completeness of the sample thanks to exten¬ 
sive simulations. Bear in mind that the simulations are run on a set of synthetic images with same 
exposure time distribution and same background of the selected SWXCS field. Another important 
aspect is that the input flux distribution of the simulated sources are taken from real, deep data. 
Point sources are ran domly extracted from a distribution modeled on the number counts found in 
deep Chandra fields ( Rosati et al. 2002b : Moretti et al. 2003 : Xue et al. 2011 : Lehmer et al. 2012 ) 
and simulated down to a flux about one order of magnitudes lower than the expected detection 
limit of the SWXCS. The flux distribution of the input extended sources was tak en fro m the num¬ 
ber counts of groups and clusters measured in the ROSAT deep cluster survey (Rosati et al. 19981) . 
Finally, to take into account the different morphologies of extended sources, we modeled the sur¬ 
face brightness of our simulated sources on real images of ten bright groups and clu sters of g alaxies 
observed with Chandra, covering a wide range in ICM temperature (for details see lLiu et al.ll2013h . 
We simulate ten times the entire SWXCS survey, which correspond to ~ 3 x 10 4 X-ray images. 


In Figure 0 we show the expected completeness (fraction of extended sources recovered at a 
given value of input net counts), and the expected number of spurious extended sources expected 
in the entire SWXCS in bins of net detected counts, in the upper and lower panel, respectively. 
We also plot the completeness and contamination separately for fields with exposure time above 
and below 50ks, which account for 10% and 90% of all the fields respectively. We note that most 
of the incompleteness and most of the spurious sources come from the 10% deepest fields. The 
reason is the less efficient performance of EXSdetect in presence of high background and crowded 
fields. For the entire SWXCS, the completeness falls below 90% at about 130 net counts, and 
reaches 83% at 80 net counts. Meanwhile, the contamination number increases rapidly below 80 
net counts. Above 80 net counts, the total number of spurious extended sources in the entire survey, 
at any flux, is estimated to be about 20, most of which with less than 150 net counts. Therefore, 
although the completeness is still high and robustly measured below 80 counts, we conservatively 
set the detection threshold of our catalog to 80 net counts within the source region as defined by 
the EXSdetect algorithm, to keep a low number of spurious detections. 


For each field, the minimum detectable count rate, computed as 80 /t exp , corresponds to a 
position-dependent flux limit obtained by multiplying this number by the field energy conversion 
factor ( ECF ) at the aimpoint in 0.5-2 keV band (which accounts for the Galactic absorption and is 
computed for and average thermal model with a temperature of 5 keV, a metal abundance of 0.3 Z G , 
and a redshift z = 0.4), and by the normalized exposure map, which accounts for the vignetting 
effects. As shown in Paper I, the ECF s depend weakly on the spectral parameters. Therefore, a 
flux-limit map is obtained for each field. The sum of the flux-limit maps of the entire set of fields 
considered in the SWXCS provides the sky coverage of the survey as a function of energy flux 
(see Paper I). The sky coverage of the SWXCS is shown in Figure 0 (solid line). The difference 
with respect to Paper I (dashed line) at low fluxes is mostly due to the lower threshold used (80 
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Measured net counts 


Fig. 3.— Upper panel: the solid red line shows the completeness of SWXCS (defined as the 
ratio of recovered sources and the number input sources) as a function of the input net counts 
as measured with our simulations. The dashed red lines indicate the completeness obtained only 
considering fields with an exposure less than 50 ks (long dashed line) and larger than 50 ks (dot 
dashed line). Lower panel: histogram of the number of spurious sources expected in the entire 
SWXCS catalog from our simulations in bins of net detected counts. The different colors show the 
number of spurious sources found in fields with exposure time above (green shaded are) and below 
(blue shaded area) 50 ks. 













- 11 - 


net counts as opposed to 100) and to the fact that several fields are now deeper thanks to new 
observations targeting fields already included in Paper I. We define the flux limit of the survey S u m 
the flux at which the sky coverage falls below 1 deg 2 . This corresponds to a flux of Su m = 7 x 10 -15 
erg cm -2 s -1 , which is fainter than that in Paper I (~ 10~ 14 erg cm* 2 s _1 ). At the bright end, the 
increase by a factor of 10 is due to the inclusion of the many shallow fields (mostly not associated 
to GRB), which were not considered in Paper I. The maximum solid angle covered by the survey 


is ~ 400 deg 2 , reached above a flux S 


3 x 10 12 erg cm 2 s l . If we compare the sky coverage of 


the SWXCS with previous deep X-ray surveys of galaxy clu sters, we find that the SWXCS reaches 
a depth similar to the ROSAT Deep Cluster Survey (Rosati et al. 19981) and a width similar to the 
ROSAT 400d survey (Burenin et al. 20071) as shown in Figure |4] (dotted and dot-dashed lines). 
Since the sky coverage is essential to derive the number counts and eventually, once the redshifts 
are available, the source number density as a function of redshift, we provide the tabulated values 
in the second column of Table Q] 



Fig. 4.— Sky coverage of the SWXCS as a function of the soft band flux (solid line). For com- 



The completeness can also be computed as a function of the energy flux, simply by computing 
the flux of each source which depends on the source counts rate, the ECF in the field and the 
and the actual effective area in the extraction region of the source (see Section 3 for details). 
The completeness as a function of the flux is then simply obtained as in the previous case, but 
computing the ratio of the recovered over the input sources in bins of energy flux. In practice, this 
is equivalent in convolving the completeness function shown in the upper panel of Figure [3] with 
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Table 1: Tabulated values of the sky coverage and of the completeness as a function of the measured 
energy flux in SWXCS. 


s 

erg cm“ 2 s 1 

Sky Coverage 
deg 2 

Ratio 

3.25e-15 

0.1 

42.5% 

5.15e-15 

0.6 

50.5% 

8.15e-15 

1.6 

60.4% 

1.29e-14 

4.4 

72.0% 

2.05e-14 

9.7 

80.5% 

3.25e-14 

18.5 

85.6% 

5.15e-14 

29.6 

90.5% 

8.15e-14 

42.9 

95.0% 

1.29e-13 

64.5 

97.6% 

2.05e-13 

104 

99.4% 

3.25e-13 

171 

100% 

8.15e-13 

319 

100% 

2.05e-12 

372 

100% 

8.15e-12 

409 

100% 
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the actual distribution of exposure time, ECF and effective area of the SWXCS as represented 
in the simulations. In Figure [5J upper panel, we compare the simulated input extended source 
distribution with the measured extended source distribution recovered with EXSdetect as a function 
of the energy flux. Note that the input distribution is given by the average of ten actual realizations 
of the input model used in the simulations, therefore it has a lcr statistical uncertainty shown by 
the shaded area. Note also that here the sky coverage is already accounted for. The curve in the 
lower panel of Figure [3] shows the ratio of output to input source distributions as a function of the 
measured flux. This function is our best estimate of the completeness of the SWXCS as a function 
of the energy flux. We find that the completeness correction is relevant below 5x10“ 14 erg cm* 2 s _1 , 
while it is negligible above 10~ 13 erg cm -2 s* 1 . We remark that, since the completeness is computed 
as a function of the measured energy flux (and not the actual input flux used in the simulation), 
this correction takes account also the effect of the Eddington bias expected in the SWXCS. This 
function is also tabulated in TableU] third column. The combination of the sky coverage and of the 
completeness function as a function of the energy flux allows one to derive the SWXCS number 
counts directly from the catalog (see Section 3). 


2.3. Filtering Of Spurious Sources Not Included In Simulations 


By running EXSdetect (version 3.0) on the soft band images of the 3000 fields we obtain 430 
extended source candidates with a soft band photometry larger than 80 net counts in the source 
region defined by the EXSdetect algorithm itself. According to the expected performance of EXS¬ 
detect, the source catalog obtained directly should have only about 20 spurious sources. Apart 
from this budget, however, we know that we have spurious sources associated to circumstances not 
included in the simulations, but that can be easily identified and filtered out. 


The first class of these sources are the piled-up targets, which we can always safely assume 
to be bright AGN or GRB. At the Swift-XRT angular resolution, in fact, it is very unlikely that 
a bright cool core can be affected by pile-up. Clearly, the pile-up effect modifies the shape of 
the inner regions of unresolved sources, causing the failure of our source classification algorithm 
based on the comparison of the source profile with the synthetic image of an unresolved source 
in the same position and with the same flux. For obvious reasons, we did not attempt to include 
this effect in our simulation. On the other hand, piled-up sources can be easily identified and 
removed. As a simple and effective criterion to identify piled-up sources, we compute the count 
rate of each source. All the data we used are observed in the Photon Counting mode, which has 
a time resolution of 2.5 seconds ([B urr ows et al. 2005 ). According to the Swift-XRT documents, 
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10 4 10 10 



Fig. 5.— Upper panel: the blue dashed line and red solid line show respectively the input number 
counts and the measured number counts from the simulation as a function of the energy flux. 
Here we used the entire set of simulations (corresponding to ten times the entire SWXCS survey). 
Shaded areas show the lcr uncertainty on the input and recovered number counts. Lower panel: 
the ratio of recovered differential number counts to the input differential number counts, with lcr 
uncertainty, as a function of the energy flux. This curve, along with the sky coverage, is used to 
compute the number counts in the SWXCS. 
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any source brighter than 0.5 counts/s in the 0.2-10 keV full band should be checked for pile-up H|. 
Taking the spectrum of 3C 273 as that of a typical unresolved source, we find that the 0.5-2 keV 
band contains ~ 50% of all photons in the full band. Thus we flag all the sources with 0.5-2 keV 
count rates > 0.25 counts/s inside the extraction region as piled-up candidates. Then, we check 
their classifications from NED, finding mostly QSOs or galaxies, while none of them, as expected, 
is associated with clusters or groups. This step allows us to reject 34 extended source candidates 
as piled-up sources. In addition, some of the GRB may escape this filter since they can suffer 
pile-up only in high flux states. Since the identification of all the GRB is straightforward, we filter 
out 18 GRBs which were mistakenly classified as extended sources. A simple visual inspection is 
performed to check whether we may find truly extended sources overlapping the GRB positions, 
but we found none. 


Nearby (z < 0.05) galaxies constitute another source of contamination for our sample. At low 
redshift, spiral galaxies appear as extended X-ray sources in Swift-XRT images, which represent 
the populations of high-mass X-ray binaries (in the case of recent starburst) or low-mass X-ray 
binaries in the galaxies. Nearby elliptical galaxies may also show X-ray emission which is related 
to hot gas in their halo. To identify them we consider the following galaxy catalogs: the Local 


(Gil de Paz et al. 

2007 

), and the third Reference Catalog of bright Galaxies ( 

de Vaucouleurs et al. 


199 lh . and select all the galaxies whose major axis diameter is larger than the half power diameter 


of the XRT PSF (18"). For these sources, the extended X-ray emission strongly overlaps with the 
optical extent of the galaxies, as seen in the Digital Sky Survey (DSS) optical images or the Sloan 
Digital Sky Survey (SDSS) images when available, and therefore is clearly dominated by stellar 
sources. In some cases, we notice some emission beyond the optical extent of the galaxies. This 
is particularly evident in galaxy pairs. However, in all the cases where the extended emission may 
be associated with hot gas around the galaxies, its emission is weak and very hard to decouple 
from the stronger X-ray emission from the disk or the bulge. In other cases, some small-extent 
diffuse emission is swamped by a central AGN. For all these sources an accurate modeling of the 
non-thermal X-ray emission is needed before the thermal component can be properly evaluated. 
Given the very small impact that such contributions would bring to our final catalog, all the 30 
extended source candidates associated with nearby galaxies are excluded. 


Another kind of spurious detections are caused by bright optical sources (generally stars) 
whose intense optical/UV emission induces significant spurious charge load in the CCD. Such 
sources are automatically screened by Swift-XRT pipeline. However, this process often leaves a 
ring-like signature which is classified as extended by our algorithm, but can be immediately spotted 


*see http://www.swift.ac.uk/analysis/xrt/pileup.php 
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by visual inspection. We find and reject 10 of these spurious detections. 

Another effect that, for obvious reasons, is not considered in the simulations, is associated to 
the wings of bright clusters which only partially fall in the FOV of the XRT image. By visually 
inspecting RASS, or XMM-Newton images when available, centered on the source position, we 
identify 11 cases in which the wings of large X-ray clusters have been detected as extended source. 
The classification of these sources by EXSdetect is indeed correct, however we will not include 
them in the catalog, since the majority of the cluster emission is, in all the cases, well beyond the 
XRT FOV. The clusters responsible of these detections are Coma, Abell 0496, Abell 1285, Abell 
1387, Abell 1767, Abell 2199, Abell 2199, Abell 2877b, Abell 3334, Abell 3395 SW, and MCXC 
J1423.8+4015. 

The X-ray images sometimes can be contaminated directly by sunlight, which creates a flare 
in the light curve and significant diffuse emission at the image borders. By checking the light 
curve of extended sources at the image borders, we find that the emission of 22 source candidates 
is actually due to optical flares. These sources are discarded as optical contaminations. 

We also find that in 11 cases our extended source candidates are associated with the position 
of the nominal target of the Swift-XRT observation within a distance of 2'. Among these, nine 
are bright ROSAT X-ray sources, which makes them possible galaxy clusters detected but not 
identified by ROSAT because of its poor spatial resolution. In other words, although the targets 
of ROSAT sources were not identified as clusters (see ^2.11) . they have higher probabilities to be 
clusters than random objects. The other two targets are galaxies which are not associated with 
known clusters. Apparently, these last two sources may be allowed in a serendipitous sample. 
However, the perfect position match (< 6") indicates a strong connection between the galaxy 
targeted by Swift-XRT and the X-ray extended source. Also in these two cases, a positive bias is 
introduced because the possibility of finding a cluster associated to the targeted galaxy. Overall, 
although these 11 targeted observations survived our field selection ( 32. lb . they may have higher 
probabilities to host groups or clusters with respect to truly serendipitous observations. Therefore 
we remove these 11 extended source candidates from the final list. 

Finally, we consider a last case which is not properly treated in our simulations. We ran¬ 
domly sampled the number counts of unresolved sources on a solid angle of 400 deg 2 . This is 
what actually happens for genuinely random fields. However, a significant fraction of fields in 
the Swift-XRT archive are targeting very bright QSO. In some of these cases we find anomalous 
extended sources which likely to be due to spurious effects associated to the X-ray emission of 
the bright QSO. In particular, in two fields targeting extremely bright quasars, we find 6 sources 
in the outskirts which are most likely associated to the anomalous background due to the presence 
of the bright source. We also find 5 very bright sources which are surrounded by a much fainter 
extended emission. In these few cases, the extended emission component is sufficient to classify 
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these sources as extended, however the central source is not identified as unresolved, probably due 
to the contribution of the extended emission, so it is not removed, as it is done in all the cases when 
unresolved sources are embedded in extended emission. These cases should be treated separately, 
and a PSF deconvolution of the unresolved source emission should be applied before the extended 
emission could disentangled from the unresolved emission and properly measured. Since this pro¬ 
cedure would introduce large errors on the photometry of the faint extended emission, we decide 
to remove these sources from the catalog, despite the fact that they do include extended emission. 
Clearly, the angular resolution is a strong limitation which hampers us to properly deal with such 
cases. Overall, we remove 11 sources due to effects associated to the presence of targeted, very 
bright unresolved sources. 

At the end of this cleaning procedure, we have removed 147 sources due to effects which 
could not be included in the simulations. In principle, these effects could be implemented in the 
reduction pipeline with some additional effort. However, in our case, due to the limited number 
of sources in the SWXCS, a manual check a posteriori is feasible and the automatization of this 
filtering process is not crucial at tis stage of the project. After this step is completed, we are left 
with 283 group and cluster candidates whose properties are well described by the completeness 
and contamination function obtained with our simulations. 


2.4. Beyond The Software: Learned Visual Inspection 


Assuming that we have filtered out any possible source of contamination not included in the 
simulations, we are now dealing with a sample with known statistical properties. However, we 
can extract more information from the simulations. If we revise all the spurious sources we find 
in the simulations (several hundreds) we can understand in most of the cases the reasons why the 
Voronoi algorithm failed. For example, the quite common case of blending of two or more visual 
sources can be easily associated to a particular pattern in the surface brightness distribution of the 
spurious source candidate. This visual learning procedure is very effective and in principle could 
be implemented in the software as a machine learning process, as is now commonly done for data 
mining in very large surveys. However, the human eye still appears to be the best tool to implement 
such complex proces ses. Th is has been clearly shown by the “ crowd-sourcing ” projects proposed 
by Galaxy Zoo (see [Li ntott e t al.1 2 0111 in the last years, which have been proved to be highly 
successful. It may be very useful, in the future, to set up such a crowd-sourcing astronomical 
project based on X-ray images. 


In our case, given the small size of our source list, we are able to train the eye with simulation 
and manually apply the filter with a direct visual inspection. The sources which are flagged as 
spurious after this step are 20, consistent with the number of spurious sources expected from the 
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simulations (~ 20). Therefore, we decide to include this step in our source classification scheme. 
In the end, we are left with a final catalog containing 263 source candidates of clusters and groups, 
whose selection function is shown in Figure [3] After the visual inspection, we are confident that 
the contamination in the SWXCS catalog is reduced at a level which can be safely ignored when 
deriving statistical properties of the sample. 


3. SWXCS CATALOG 



Fig. 6.— Upper panel: the red solid line shows the corrected number counts of SWXCS. The 
black line is the best-fit double power-law model of RDCS ( Rosati et alll998 ). Shaded areas show 
the corresponding lcr confidence intervals in both panels. Lower panel: the ratio of SWXCS over 
RDCS differential number counts. 
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In Table [2] we list the 263 sources of the catalog with their X-ray properties. The sources 
already presented in Paper I are marked with an asterisk. The catalog contains the following 
information: 


• Column 1: source name according to the format officially accepted by the IAU Registry in 
February 2013. The format is SWXCS JHHMMSS+DDMM.m. This format is different 
from that used in Paper I, but it was already used in Paper II. Note that for the sources that 
presented in Paper I, we keep the same positions used in Paper II, although the new positions 
typically differ by ~ 5". 

• Columns 2 and 3: RA and DEC coordinates of the X-ray centroid, defined as weighted 
median position of the 27 brightest pixels in the source region (each pixel is weighted by its 
density which equals to the pixel value divided by Voronoi cell area). 


• Column 4: the effective exposure time of each source computed as: 

, _ Z n i 

*eff ~ v n It ’ 

Zj fti! 


(i) 


where i is the index of the filled pixels within the source extraction region, «, is the photon 
count in the i th pixel, U is the corresponding value in the vignetted exposure map. 


• Column 5: value of the Galactic neutral hydrogen column density in unit of 10 
found in the Leiden/Argentine/Bonn radio survey (Kalberla et al. 20051) . 


,20 cm 2 , 


as 


• Column 6: R e ff in arcsec, defined such that /r/C /7 . equals the source extraction area (which 
has no a priori constraints on its shape). 


• Column 7: the net counts N net measured in the soft band image by EXSdetect inside the 
source extraction region found by the VT+FOF method, after background subtraction and 
removal of unresolved sources. At present, no correction is applied to compensate for the lost 
diffuse emission in the region of the removed unresolved source. The lcr error is computed 
as -^N t ot + Nbkg, where N hkg is the counts of background photons, and N tol = N net + N hkg . 

• Column 8: signal-to-noise ratio (S/N) in the soft band, computed as the net counts divided 
by the associated error. 

• Column 9: estimated soft band flux in units of 10“ 14 erg cm" 2 s" 1 . To compute the flux, 
for each field, we estimate the energy conversion factor ( ECF ) in 0.5-2 keV band, taking 
into account the Galactic absorption in this field, assuming a hot diffuse gas emission model 
with a temperature of 5 keV, a metal abundance of 0.3 Z Q , and a redshift z = 0.4. As shown 
in Paper I, the ECF depends weakly on the spectral parameters. The flux of each source 
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is calculated as 5 = N net x ECF/t e ff. The lcr error of flux is measured by propagation 
considering the error of net counts and a 4% systematic error of ECF due to the typical 
uncertainty in the actual spectra shape of each source (see Paper I). More accurate fluxes 
will be available from the X-ray spectral analysis of our sources, which is postponed to a 
forthcoming paper (Moretti et al. in preparation). 


Finally, using the sky coverage and the completeness function described in ^12.21 we compute 
the number counts of the SWXCS catalog, following the same procedure as described in Paper 
I. The corrected number counts are shown in the upper panel of Figure [6} with lcr confidence 
intervals, which includes the Poissonian error and the uncertainties on the average conversion 
factors. As we also found in Paper I, the_numbercounts are consistent with the logN-logS of the 
ROSAT Deep Cluster Survey (RDCS) ( Rosati et al. 1998 ) (see lower panel of Figure [6]). The faint 
end is also consistent with the very-deep number counts measured from the Extended Chandra 
Deep Field South (Finoguenov et al. in preparation). We find that the differential number counts, 
after correction for incompleteness, can be fit with a broken power-law. The best-fit model is: 


dN_ 

~dS 


(1.29 ±0.26) x 


(1.75 +0.35) x 


- 2 . 14 ± 0.06 


10 - 13 , 

£ x - 1 . 68 ± 0.16 


10 


-13 


S > (5.1 +2.0) x 10 


S < (5.1 +2.0) x 10 


-14 


-14 


( 2 ) 


where S in unit of erg cm 2 s 1 . Therefore, the slope of the faint end appears to be flatter than the 
slope at the bright end, although with a low significance^ 3cr). 


4. COMPARISON WITH PAPER I 

The first release of the SWXCS catalog (Paper I) was based on a much smaller number of 
fields (~ 300 GRB follow-up fields), and was obtained using a standard wavelet detection algo¬ 
rithm coupled with a growth-curve method used to characterize extended sources. To compare 
this work with the previous release, first we investigate the effect of using EXSdetect, which is 
applied to real data for the first time in this work. After filtering out all the spurious sources as 
described in ^2.31 and ^2.41 EXSdetect detects 113 sources in the GRB follow-up fields which 
were used in Paper I (and clearly also included in this work). All of the 72 sources presented in 
Paper I except one (SWXCS J022344+3823.2) are recovered by EXSdetect. In addition, 42 new 
sources are detected for the first time by EXSdetect. Among the newly detected sources, 17 have 
less than 100 net counts. The other 25 new sources, instead, have a photometry brighter than 100 
net counts, and therefore should have been included in the first release of the catalog in Paper I. 
So we conclude that they were simply missed by the detection method used in Paper I. This shows 
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that the EXSdetect algorithm is more efficient, allowing us to recover, above the same photometry 
threshold, a number of sources 30% higher than in Paper I. In the upper panel of Figure [71 we plot 
the flux distribution of the 113 EXSdetect detected sources, compared with the flux distribution 
of the sources of Paper I. Most of the newly detected sources are found at fluxes below 10" 13 erg 
cm 2 s _1 , showing that EXSdetect is able to reach higher sensitivity allowing us to further explore 
the flux range where medium and high-z clusters are found ( few x 10" 14 erg cm" 2 s" 1 ). 

In the lower panel of Figure [7] we show that the photometry obtained with EXSdetect is in 
good agreement with the values found in Paper I for the 71 sources in common. Note that the 
fluxes measured in Paper I are corrected for the missed flux beyond R ext assuming a best-fit [3 
model, while in this work no correction is applied. The best-fit relation between EXSdetect fluxes 
and the fluxes in Paper I reads: 

l0g (l^) = (0-97 ± 0.13) xlog(^§)-0.05. (3) 


The EXSdetect fluxes are somewhat higher than the fluxes measured in Paper I, up to a maximum 
of ~ 10% for the brightest sources. On the other hand, our simulatio ns sh owed that the EXSdetect 
photometry is accurate at the level of ~ 1 - 2% (see Figure 10 in Liu et al. [20131) . We remind 
the reader that EXSdetect automatically defined irregular extraction regions thanks to the Voronoi 
algorithm, as opposed to the circular extraction regions defined in Paper I. We find that generally 
the EXSdetect extraction region is larger than the circular region used in Paper one, as shown in 
Figure [8j by comparing the R e ff with the extraction radius R ext which was defined in Paper I as 
the radius where the average source flux equals the background level. Note that 2 sources whose 
R e ff is more than twice larger than R ext do not appear in this figure, because of large, low surface 
brightness extents associated to these sources, which was not accounted for in Paper I. Therefore, 
we conclude that the use of extraction region defined by the Voronoi method is more efficient in 
recovering the flux in the low surface brightness outs ki rts of extended sources, providing a more 
accurate estimate of the total flux. 


Finally, for the sources in common, we show the distribution of the displacements between 
the positions published in Paper I and the positions measured by EXSdetect (see Figure [9]). The 
discrepancy is mostly due to the different definition of the center of the source used in this work 
and, to a lesser extent, to the larger extraction regions. Despite this, the center of the large majority 
of the sources is changed by less than 20", a value very close to the HEW of the Swift-XRT PSF. 
Eight sources are found at separations between 0.5'and 2', because of the large extent and the rather 
flat surface brightness distribution of these sources. Although the change in the nominal position 
of some source would imply a change in the name according to the IAU format, we prefer to keep 
the same name used in Paper II for the sources of the first release of the SWXCS catalog. 
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S erg cm 2 s 1 , EXSdetect 


Fig. 7.— Upper panel: histogram distributions of the EXSdetect measured soft-band fluxes 
of all the sources detected with EXSdetect in the same GRB fields as used in Paper I (red 
histogram) and of the 71 SWXCS sources that already presented in Paper I (blue filled his¬ 
togram). Lower panel: for the 71 SWXCS sources that Paper I has in common, the soft-band 
fluxes S p aper! measured in Paper I within R ext , compared with the soft-band fluxes S E xs mea¬ 
sured in this paper with EXSdetect from a source region defined by the VT+FOF algorithm. 
Black solid line shows the relation S p ap eri = S exs , while the blue dashed line shows the best 
fit log(Sp^.,/10- 13 ) = 0.97 x logGS£* 5 /FT 13 ) - 0.05. 
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Fig. 8.— Extraction radius R ext used in Paper I to define the source region, compared to the 
effective radius R e ff for EXSdetect sources. The solid line shows the relation R ext = R e ff. 


E 



Fig. 9.— Displacements between the source positions in Paper 1 and the positions found in this 
work, in units of arcsec. The green circle shows the median separation corresponding to 5". 
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5. CROSS-CORRELATION WITH OPTICAL, X-RAY AND SZ CATALOGS AND 

OPTICAL LOLLOW-UP 


We checked for counterparts in previous X-ray cluster surveys, in optical cluster surveys, and 
in the Planck SZ cluster survey. We simply assume a search radius of 2' from the X-ray centroid, 
which has been shown to be an efficient criterion in Paper I. Nevertheless, we also inspected the 
area within 5' from the X-ray centroid to investigate whether some possible identification is found 
at radii larger than 2'. Counterparts at distance between 2' and 5' are included when the optical 
or SZ corresponding source has a large uncertainty in position. This is often the case for optical, 
sparse clusters, or for SZ cluster candidates. We list all the counterparts associated to the SWXCS 
sources in Table [3j with the measured redshift when available. In case of multiple counterparts, 
we list all of them. Except for a few cases where we have multiple counterparts with statistically 
inconsistent redshifts, we keep the counterpart with the smallest distance from the X-ray center. 


From optical surveys, we found 233 optical counterparts corresponding to 116 SWXCS sources, 
including 84 from the SDSS WHL catalog ( Wen et all 2012). 25 from the SPSS AMF catalog 
( Szabo et al. 2011), 28 from the S PSS MaxBCG catalog ( Koester et al. 2007), 45 from the SDSS 
GMB CG catalog ( Hao et al.ll2010 ) , 8 from the SDSS C4 catalog ( Miller et al. 12005 : von der Linden et al 


20071) . 27 from the Abell c atalog (lAbell et al.lll989l) . 8 from the NSC S1 catalog (IGal et al 


4 from the NSCS2 catalog (Gal et al. 20031) . 3 from the EDC catalog ( Lumsden e t al. 


at last, 1 from the SDSS galaxy groups and clusters catalog built by IBerlind et al.! ( 2006 ). The 


2003D, 


19921). and, 


majority of the SWXCS sources with optical counterparts are listed in more than one catalogs. A 
few WHL counterparts published in Paper II are found with different names in this work, because 
of the updated version of WHL catalog used here. 

From X-ray surveys, we found 70 X-ray counterparts classified as cluster, corresponding to 36 


SWX CS sources. In detail, we found 12 X-ray clusters in the RO SAT 400d catalog (IBurenin et. al. 


2007b, 11 in the Northern ROSAT All-Sky (NORAS) catalog (Bo hringer et al.112000 ). 8 in the 


ROSAT-ESO flux Limited X-ray galaxy cluster c atalog (REFLEX iBohringer et al. 


the XMM-Newton Cluster Survey (XCS) catalog (Mehrt ens et al 


2012), and 1 in the Chandra 


», 


3 in 


Multiwavelength Project (ChaMP) gala xy cluster catalog (IBarkhouse et al.1120061) . We also found 
35 counterparts in the MCXC catalog (Pi ffarctti et al . 2 0111) . which includes most of the X-ray 
clusters above. 


Finally, for 15 SWXCS sources, we foun d 16 cluster counterparts detected v ia SZ effect 


13 by Planck (IPlanck Collaboration et al.ll2013f) and 3 by South Pole Telescope (SPT. lBleem et al. 


2014). The Planck sources are typically at larger distances from the X-ray centroid (bet ween Y and 


3'), b ecause of the much larger position errors of Planck clusters (see iMalte Schafer & Bartelmann 


20071) . 
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Overall, about half (137) of the 263 SWXCS sources have been previously identified as groups 
or clusters of galaxies, while 126 SWXCS sources are new cluster and group candidates. Thanks 
to these identifications, we are able to recover the redshift information for a significant fraction of 
our sample. We collect spectroscopic or photometric redshift for 130 of our sources. Moreover, 
to increase the number of available redshifts, we also search in NED catalogs for single galaxies 
with published redshift not associated to previously known clusters within a search radius of 1" 
from the X-ray centroid of our sources. We find 50 galaxies with measured redshift for 47 of our 
sources, as a complement to the redshifts obtained from cluster counterparts. In 35 cases where we 
have both cluster and galaxy counterparts, the galaxy redshifts are consistent with those of clusters. 
In the 12 cases where no cluster counterpart is found, we tentatively assign the galaxy redshift to 
our X-ray source. 

If we consider also the X-ray redshift derived from the Telescopio Nazionale Galileo (TNG) 
observations and the X-ray spectral analysis of the sources in catalog I (Paper II) we have a total 
of 158 sources with redshift, from optical spectroscopy or photometry, or from X-ray spectral 
analysis. Therefore, about 60% of our sample has redshift information. For these sources, we plot 
the rest-frame 0.5-2 keV luminosities versus redshifts in Figure [TOl 



z 


Fig. 10.— Rest-frame 0.5-2 keV luminosity versus redshift. More accurate luminosities measured 
with X-ray spectroscopy in Paper II are used when available (red points). 

We remark that 116 sources overlap with SDSS images. In FigureQT]we show, as a sample, 
a selection of SDSS r-band images of SWXCS sources with obvious optical counterparts, with 
X-ray contours overlaid. The X-ray and SDSS images (when available) for all the SWXCS sample 
can be found in the websites SWXCS website (http://www.arcetri.astro.it/SWXCS/ and 
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http: //swxcs. ustc. edu. cn ). Among the source with optical redshift, the highest redshift is 
z = 0.92 for XMMXCS J142908.4+424128.9 ( Mehrtens et al. 2012 ). This confirms that the depth 
of our catalog is sufficient to select clusters up to z ~ 1. Another high redshift source candidate 
is SWXCSJ011432-4828.4, whose redshift is measured to be 0.97 + 0.02 from X-ray spectral 
analysis in Paper II. The presence of clusters at z ~ 1 is not unexpected in the SWXCS, given 
the non negligible sky-coverage of the SWXCS at low fluxes. Based on previous results from 
the RDCS, we expect of the order of ~ 10 clusters with z > 1. In addition to the few sources 
already mentioned, we already identified a sample of high-redshift candidates among the sources 
with SDSS images but no optical counterparts. As an example, the SDSS images of four of our 
high-z cluster candidates are shown in Figure [T2l 



SV?X£S J152252*2535.6 z=0.5S SWXCS JO92650.3013.8 z=0.56 
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SWXCS J132605+1223.0 z=0.21 SWXCS JJ 52343+5951.8 z=0.36 


SWXCS J140637.2743.6 x=0.65 SWXCS ^063340.3311.0 1=0.71 


Fig. 11.— SDSS r-band images of a few medium-z SWXCS sources with obvious optical coun¬ 
terparts. The images have sizes of 5' x 5' or 10' x 10'. The X-ray contours (green lines) correspond 
to 2, 5, 10, 30, 70, 150, 300 times the local background. 


The next step of our project is to increase the number of identifications and redshift measure¬ 
ments, in order to use our sample for statistical studies and cosmological tests. We have started an 
extensive follow-up program to obtain sensitive, multi-band imaging photometry of the SWXCS 
sample. Our immediate goal is to measure the integrated properties of the stellar populations 
of the galaxies through SED fitting. We will also explore the correlation of the galaxy proper¬ 
ties with those of the hosting cluster and their evolution with redshift. The planned observations 
consist of deep CCD images in the UBVRIz filters which, in the redshift range 0.3 < z < 1, 
probe the rest-frame spectral range from mid- and far-UV to the optical wavelengths. This choice 
allows us to properly sample the wavelengths across the 4000 A break, which are key for accu¬ 
rate measurement of the integrated stellar mass, star-formation rate, average dust obscuration and 
luminosity-weighted age through SED fitting. For the southern sources we use the Du Pont 2.5- 
meter telescope at Las Campanas Observatory coupled with the Cassegrain-focus Direct Camera. 
Sources in the northern emisphere are observed using the Mayall 4-meter telescope with the MO- 
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SWXCS J143131+3622.9 ! SWXCS J143223+3617.9 


Fig. 12.— SDSS r-band images of four SWXCS sources with soft-band fluxes less than 5 x 10 -14 
erg cm' 2 s^and without obvious optical counterparts. These sources are among the high-z cluster 
candidates in the SWXCS. The image sizes are 5' x 5'. The X-ray contour generating method is 
the same as used in Figure [TTl 
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SAIC Prime Focus camera at the Kitt Peak National Observatory. Currently, we have observed a 
total of 41 groups and clusters, 11 in the South and 30 in the North one. We plan to release the 
reduced and calibrated images and source catalogs of the first year of observations in early 2015. 
The program will continue in the following years. 


We finally note that we can add a significant number of redshifts by extending the X-ray 
spectral analysis to the entire sample. Although t he req uirements for a successful identification of 
the redshifted K a Fe line, as shown in Yu et al . (20UJ) for Chandra, do not apply to most of the 
SXCS sources, the lower background and the slightly better spectral resolution of XRT allows X- 
ray redshift measurements in a lower S/N regime, as shown in Paper II. The X-ray spectral analysis 
of the SWXCS sample will be presented in a forthcoming paper (Moretti et al. in preparation). 


6. CONCLUSIONS 

We search for candidates groups and clusters of galaxies in ~ 3000 extragalactic Swift-XRT 
fields. These fields are selected in order to provide a truly serendipitous survey, therefore excluding 
all the fields in the Swift-XRT archive which are somehow correlated with galaxy clusters and 
groups. We use the software EXSdetect, which has been specifically developed for this project and 
it is optimized for detection and photometry of extended sources in Swift-XRT images. Therefore, 
both in terms of covered solid angle and of sensitivity, this work constitutes a significant extension 
of the first SWXCS catalog published in Paper I. 

We find 263 X-ray extended sources (including the 71 sources already presented in Paper I) 
with negligible contamination and a well defined selection function. The sky coverage ranges from 
a maximum of ~ 400 deg 2 to 1 deg 2 at a flux of 0.7 x 10~ 14 erg cm“ 2 s -1 . The logN-logS is in very 
good agreement with previous deep surveys. We cross-correlate SWXCS sources with previously 
published optical, X-ray or SZ cluster catalogs, finding that 137 sources are already classified as 
clusters in any of the three bands, while 126 sources are new cluster and group candidates. We al¬ 
ready collected redshifts measurement (both optical, spectroscopic of photometric, and X-ray) for 
158 sources (60% of the sample). When the optical follow-up and the extension of the X-ray spec¬ 
tral analysis will be completed, the SWXCS will provide a large and well defined catalog of groups 
and clusters of galaxies to perform statistical studies of cluster properties and tests of cosmological 
models. All the results of the SWXCS are publicly available on http://www.arcetri.astro.it/SWXCS 
or http://swxcs.ustc.edu.cn, including machine-readable tables and the EXSdetect code. 
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APPENDIX 


A. Background Estimation 


The background of an X-ray image is defined as the sum of all the recorded photons not 
associated to astronomical sources, or associated to some astronomical component that can not be 
resolved (like the Galactic diffuse emission). Practically, we divide the photons in an X-ray image 
into two components: a background component with a roughly constant flux distributed randomly 
across the whole field, and an additional component, associated to single sources, with a highly 
concentrated spatial distribution covering only a very minor fraction of the field. A well known 
result, obtained by numerical simulation (Kiang 1966 ) is that for randomly positioned points, the 
distribution of the Voronoi cell areas follows an empirical formula: 


P(f) 


A/f 32 8 4 

4// (—- + — + — + 1 ), 

3/ 3 f 2 f 


(A.l) 


where P is the cumulative probability distribution function (CDF, P e [0,1]), / = 1/a is the 
inverse of the cell area a, called flux here, and / = //(/) is the flux normalized to the average value 
(/) = 1 /(a). Therefore, the distribution of the Voronoi cell areas is provided by a function of / with 
only one parameter (/). The value of (/) is the background of the image. This relation was used in 
the first Voronoi algorithm for X-ray source detection proposed by Ebeling & Wicdcnmann (19931) . 
Photons in the faintest end of the filled-pixel distribution can be assumed to be only due to the 
background. Therefore, they proposed an accurate measurement of the average background flux 
by fitting the faint-end CDF, specifically, only for the //(/) < 0.8 part. This part includes ~ 27% 
of all the filled pixel in a pure-background image. We made use of this method in EXSdetect in 
Liu et al. (120130 . Here we improved the background estimation running Monte Carlo simulations 
as described below. 


We start from a direct test of Equation IA.II with simulations. We randomly distribute one 
million photons in images of different sizes, chosen in order to have three different flux levels 
(0.001,0.03,0.1 photon/pixel). Note that float-value positions are assigned to each photon, which 
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is equivalent to assuming an infinitely small pixel scale for the images. Using the SweepLine 
subtask of EXSdetect, we construct the Voronoi diagram for each image and calculate the area 
of each cell. Then we compare Equation IA. II with the simulated distribution of fluxes. We find 
that Equation IA. 1 1 describes the CDF well in the entire range, but not particularly in the faint end. 
Therefore, the best fit of Equation IA.1I for the //</) < 0.8 part does not recover the average 
background flux accurately. A systematic deviation is introduced. 

Further more, the assumption of infinitely small pixel scale is not realistic. Photon positions 
in real images are always in integer rather than in float. In other words, digital images always 
have limited resolutions. The limited resolution induces two effects: a distortion of the CDF in the 
faint-end which is due to the lower limit of Voronoi cell area, and a reduction of data points which 
is due to the fact that multiple photons on the same pixel contribute a single Voronoi cell. These 
effects are flux dependent, being more significant at high fluxes. In the simulation, we convert the 
photon positions into pixel positions, and test Equation lA.ll in this realistic situation. As shown 
in Table IA.fl the usage of Equation IA. 1 1 underestimate the average flux by 3-6% in the flux range 
0.001-0.1 photons/pixel. 

The background measurement is improved as follows. We add a parameter c into Equation 

EH 


P(f,c) = e-'rt-C + S- + £ + 1). (A.2) 

2 /- / 

Equation IA.2I corresponds to Equation IA. 1 1 when c - 4. Then we repeat the fit of the simu¬ 
lated CDF in the region //</) < 0.8 and search for value of the c parameter which minimize the 
difference between the model and the simulated data. We also replace the definition of average 
flux (f) = 1 /(a) = N P i xe i s / X a with (/) = Nphotons / X a ' to ta ke into account multiple photons in 
the same pixel. As shown in Figure lATl we find a linear correlation between c and (/): 


c = -0.63(f) +3.94. 


(A.3) 


Equation IA.2I and IA.3I allow us to recover the background accurately with no deviation, as 
shown in Table lA.il 

A main limitation of our background measurement method consists in the background vari¬ 


ation across the image due to vignetting. In order to take vignetting into account, in Liu et al. 


(i2013h . the background map of the field, in units of photon/pixel, was obtained by multiplying 
the average background flux by the vignetted exposure map. However, only the background from 
astronomical sources is vignetted, the background components associated to instrumental noise 
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Table A.l: Background estimation accuracy, measured by the ratio (percentage) of output flux to 
input flux in simulation. Three levels of background flux (in unit of photon / pixel) and three cases 
of exposure depth (different number of photons within one image) are considered. 


Flux 

N=100 

N=1000 

N= 100000 


Kiang (1966) formula 

0.001 

96.5+6.3 

97.4+2.2 

97.4+0.2 

0.03 

95.9+6.4 

96.8+2.2 

96.9+0.2 

0.1 

94.2+6.1 

95.5+2.0 

95.7+0.2 


improved formulae 

0.001 

98.6+6.8 

99.7+2.3 

99.8+0.2 

0.03 

98.7+7.0 

99.8+2.3 

100.0+0.2 

0.1 

98.6+6.8 

100.0+2.2 

100.1+0.2 



Fig. A.l.— The correlation between the best fit c parameter in Equation IA.21 and the input back¬ 
ground. The red line is the best fit Equation IA.31 
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and cosmic ra ys v ary following a more constant pattern. Therefore, the simple procedure adopted 


in 


Liu et al. (120131) over-corrected the vignetting effect, especially at the image borders. 


In the new version of EXSdetect used in this work, we refined the background estimation 
in the following way. First we divide the field into about five concentric regions delimited by 
the smoothed contours on the exposure map (see upper panels of Figure lATZl) . We calculate the 
average background flux in each of these regions with the method based on the improved iKiang 
( 19661) formula, creating a step-like background map. Using the background fluxes and exposure 
times in these bins, we interpolate the relation between the background and the exposure time 
values in each concentric region with a linear regression. Applying this relation to the original 
exposure map directly provides a continuous background map (lower panels of Figure IA.2I) . The 
background map is used to recover the background value at the source position when applying the 
criterion to define source region. 



Fig. A.2.— Top left: exposure map of a SWXCS field; top right: concentric regions are obtained 
according to the exposure map values; bottom left: the background is computed in each concentric 
region by applying Equation lA.21 bottom right: a continuous background image is obtained. 


B. Source Classification 

High angular resolution is a crucial parameter for detection of extended sources in X-ray 
images. With HEW of the order of 5", it is possible, in principle, to identify any extended source 
associated to clusters and groups up to high redshift. The only critical aspect left in high resolution 
X-ray images concerns the discrimination of extended features with very low surface brightness 
from background fluctuations. In Wolter Type I X-ray mirrors, the angular resolution is maximized 
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at the aimpoint, while at large off-axis angle the PSF rapidly degrades. This aspect creates several 
problems when searching for serendipitous extended sources across an image, particularly if the 
PSF is not axisymmetric and if the image is obtained by the merging of many exposures. 


Swift-XRT has the valuable property, unique among existing X-ray facilities, of showing a 
constant PSF across the FOV, but the price to pay is a moderate angular resolution. In many cases, 
the extent of the image of an unresolved source is not very far from the extent of a genuine, compact 
extended source, like clusters at high-z, cool core small clusters, etc. To minimize this effect, we 
significantly improved the source classification method of EXSdetect with respect to the version 
of Liu et al. ( 2013 ). 


In Liu et al. (2013 1 the source image and the PSF model were compared inside a circle with a 
radius of 5 pixels (corresponding to a ~ 60% encircled energy at 1.5 keV) to take advantage of the 
larger S/N ratio with respect to the outer part of the PSF. In other words, the difference in the profile 
of the extended source image with respect to the PSF model within 5 "was sufficient to identify it as 
extended. Sources were then divided into three types: I unresolved sources; II ambiguous sources; 
III extended sources. These selection thresholds were obtained based on our simulations, as shown 
in Figure IB.ll In this Figure we show the value 1 — P where P is the probability of a source of 
being extended, corresponding to our selection threshold as a function of the S/N. Color coded is 
the density, in the 1 - P-S /N space, of the simulated unresolved (upper panel of Figure IB.ll) and 
of the simulated extended (lower panel of Figure IBTI) sources. 


In the most recent version of EXSdetect, we introduce a further step to classify the sources 
lying between the unresolved and extended regions. We argue that, in addition to the profile of the 
surface brightness distribution, an additional information is contained in the shape of a source. In 
particular, any unresolved source is expected to have approximately circularly symmetric isocon¬ 
tours, according to the PSF model. To recover this information we consider a larger radius of 7 
pixels (corresponding to ~ 70% encircled energy for an unresolved source). The PSF model clearly 
provides the flux level of an unresolved source at a radius of 7 pixels. Then, a FOF algorithm is 
run on all the pixels whose flux is larger than the threshold. If the source is truly unresolved, this 
region should be very close to a circle with a radius of 7 pixels. Significant emission detected 
outside this radius, is taken as a hint of an extended source. On the other hand, if several pixels fall 
below this value within the circle of 7 pixels, the source is most likely unresolved. 

A few cases of source disambiguation are shown in Figure IB.2I Unresolved sources can be 
identified with this method even under the contamination of other nearby sources, both brighter and 
fainter (upper panels of Figure lR2l) . On the other hand, extended sources may by misclassified 
for several reasons: some may show a very low S/N in the core; some may harbor very compact 
cores; some others simply harbor a bright unresolved source embedded in the diffuse emission. 
All these cases can be identified simply applying our disambiguation criterion (see lower panels 
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of Figure IB.2I) . We tested against simulations that this criterion is efficient when the extended 
emission is above the background across a region with an effective radius R e ff ~ 33". Clearly, the 
angular resolution of the instrument constitute a hard limit below which extended sources can not 
be identified by any mean. In the case of SWXCS, all the extended sources with R e ff close to the 
hard limit set by the HEW = 18"can not be identified as extended. Indeed, the minimum size of 
the sources in the SWXCS corresponds to R e f f = 27". 



Fig. B.l.— Source classification curves in the 1 - P-S/N space, where P is the probability of a 
given source to be extended, based on a comp arison between the surface brightness distribution 
and a PSF model within 5 pixels (see Liu et al. (120130 for details). The selection criteria (blue and 
red fines) classify the sources into three categories: I = unresolved; 11= ambiguous; 111= extended 
sources. The color-coded grid shows the number of simulated source in the 1 - P- S/N space, 
separately for unresolved (upper panel) and extended (lower panel) sources. 
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Fig. B.2.— Four examples of disambiguations in real XRT images. Red circles mark the position 
of the candidate extended source, and have a radius of 7 pixels. Green lines enclose the FOF 
regions which reach the same flux level as the PSF at the radius of 7 pixels. Blue circles mark 
pixels outside the FOF regions but inside the red circles. The two sources in the upper panels are 
classified as unresolved, while the sources in the lower panels are classified as extended. 
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Table 2. SWXCS catalog. Sources marked with an asterisk are included in the first relaease of 
SWXCS (Paper I; Paper II). Column 1: source name according to the IAU Registry; sources 
included in the first release keep the name used in Paper II despite the new centroid positions. 

Column 2-3: RA and Dec (J2000) coordinates of the X-ray centroid in degree. Column 4: 
effective exposure time at the source position in sec. Column 5: Galactic HI columns density in 
10 20 cm -2 . Column 6: effective radius R e /f in arcsec; the area of the source region is /rR 2 /r 
Column 7: net counts in the source region in the 0.5-2 keV band with 1 cr error. Column 8: S/N in 
the 0.5-2 keV band. Column 9: flux in units of 10 -14 erg cnD 2 s~'in the 0.5-2 keV band, with 1 cr 

error. 


Name 

Ra 

Dec 

teff 

N h 

R eff 

Nnet 

SNR 

Flux 

SWXCS J000251-5258.5 

0.713858 

-52.974476 

304043 

1.59 

39.3 

104+17 

6.3 

0.8+0.1 

SWXCS J000315-5255.2* 

0.813067 

-52.915205 

313159 

1.59 

79.8 

1089+43 

25.5 

8.5+0.5 

SWXCS J000324-5253.8* 

0.846572 

-52.899251 

294650 

1.59 

90.5 

1012+43 

23.3 

8.4+0.5 

SWXCS J000345-5301.8* 

0.934069 

-53.030934 

308065 

1.61 

67.5 

362+30 

12.2 

2.9+0.3 

SWXCS J002044-2544.0 

5.181310 

-25.733882 

3272 

2.43 

217.2 

237+19 

12.5 

180.8+16.2 

SWXCS J002114+2059.7 

5.309518 

20.995620 

95736 

3.86 

69.9 

90+14 

6.4 

2.4+0.4 

SWXCS J002437-5803.9* 

6.157687 

-58.064728 

71839 

1.22 

114.4 

360+24 

14.7 

12.2+1.0 

SWXCS J002729-2326.1 

6.870954 

-23.435391 

28798 

1.75 

99.2 

104+13 

7.9 

8.9+1.2 

SWXCS J002824+0927.1 

7.098729 

9.451710 

44478 

3.95 

99.8 

121+15 

7.9 

7.1+0.9 

SWXCS J002826+0918.3 

7.110244 

9.304376 

45409 

3.76 

102.4 

140+16 

8.6 

8.0+1.0 

SWXCS J003316+1939.4* 

8.319273 

19.656753 

44727 

3.98 

76.1 

93+13 

7.3 

5.4+0.8 

SWXCS J003759-2504.4 

9.494119 

-25.073476 

7612 

1.45 

217.7 

94+14 

6.8 

30.1+4.6 

SWXCS J004311-1129.3 

10.796565 

-11.488947 

9586 

2.38 

100.2 

87+11 

8.3 

22.7+2.9 

SWXCS J005059-0929.5 

12.744153 

-9.491900 

30035 

3.51 

131.6 

255+20 

12.8 

21.9+1.9 

SWXCS J005500-3852.4* 

13.750356 

-38.874722 

40510 

3.31 

106.2 

136+17 

8.1 

8.6+1.1 

SWXCS J010030-4749.4 

15.126057 

-47.823462 

15702 

1.89 

129.2 

112+14 

8.2 

17.7+2.3 

SWXCS J010955-4555.9 

17.481203 

-45.930912 

198459 

1.92 

484.5 

40091+257 

156.1 

498.5+20.2 


Note. — This table is available in its entirety in machine-readable form. 






Table 3. Catalog cross-correlation results. (1) Source name; (2) optical redshift of cluster or galaxy counterparts, (p) for 
photometric; (3) TNG measured redshift from Paper II; (4) X-ray redshift from Paper II; (5) catalog where the cluster 
counterpart is from; (6) cluster counterparts; (7) galaxy counterparts within 7 arcsec found in NED; (8) separation of the 

matches in arcmin. 


Name 

Zopt ZTNG 

zx 

Catalog 

Cluster counterpart 

Galaxy counterpart 

Separation 

SWXCS J000315-5255.2* 


0.62 ±0.1 




0.00 

SWXCS J000324-5253.8* 


0.76 ± 0.01 




0.00 

SWXCS J002044-2544.0 

0.141 


Planck 

Planck 119 


1.22 


0.1410 


REFLEX 

RXC J0020.7-2542 


1.45 


0.1410 


MCXC 

MCXC J0020.7-2542 


1.46 


0.142352 


Abell 

Abell 0022 


2.08 

SWXCS J002437-5803.9* 


0.195 ±0.012 




0.00 

SWXCS J002824+0927.1 

0.2238(p) 


WHL 

WHL J002827.3±092612 


1.26 

SWXCS J002826+0918.3 

0.2258(p) 


WHL 

WHL J002826.9+091824 


0.17 

SWXCS J003759-2504.4 

0.063600 


Abell 

Abell 2800 


0.68 

SWXCS J005059-0929.5 

0.2000 


WHL 

WHL J005058.1-092923 


0.19 


0.19997 


GMBCG 

GMBCG J012.74197-09.48965 


0.19 


0.200058 


MaxBCG 

MaxBCG J012.74197-09.48965 


0.19 


0.199 


400d 

400d J0050-0929 


0.33 


0.1990 


MCXC 

MCXC J0050.9-0929 


0.34 

SWXCS J005500-3852.4* 



EDCC 

EDCC 493 


0.96 


0.164127 


NED 


LCRS B005239.6-390844 

0.11 

SWXCS J010955-4555.9 

0.0238 


MCXC 

MCXC J0110.0-4555 


0.84 


0.0238 


REFLEX 

RXC J0110.0-4555 


0.87 


0.024700 


Abell 

Abell 2877 


2.44 


0.0238 


Planck 

Planck 1024 


3.29 

SWXCS J011432-4828.4* 


0.97 ± 0.02 




0.00 

SWXCS J012457-8104.9 

0.063377(p) 


Abell 

Abell S0158 


0.13 


0.087637 


NED 


2MASX J01245648-8104579 

0.03 


Note. — This table is available in its entirety in machine-readable form. 

















